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We study the impact of the nuclear dependence of R = ctl/o"t on the extraction of the F A /F^ 
and F A /F® structure function ratios from the data on the a A ja D cross section ratios. Guided 
by indications of the nuclear dependence of R from the world data, we examine selected sets of 
EMC, BCDMS, NMC and SLAC data and find that F^/Ff < a A /a D < Ff/Fj 3 . In particular, 
we observe that the nuclear enhancement (antishadowing) for F A /F® in the interval 0.1 < x < 0.3 
becomes significantly reduced or even disappears, which indicates that antishadowing is dominated 
by the longitudinal structure function Fl- We also argue that precise measurements of nuclear 
modifications of R and F A have the potential to constrain the poorly known gluon distribution in 
nuclei over a wide range of x. 



INTRODUCTION 



Since the early lepton scattering experiments that dis- 
covered the substructure of the nucleon and eventually 
led to the development of Quantum Chromodynamics 
(QCD) as the theory of the strong interaction, deep in- 
elastic scattering (DIS) has been a critical tool in the 
investigation of the quark and gluon structure of nucle- 
ons and nuclei. While initially nuclear effects in DIS 
were thought to be largely negligible, this was proven 
wrong by the measurement of the ratio of the iron to 
deuterium structure functions performed by the Euro- 
pean Muon Collaboration (EMC) at CERN in 1983 Q. 
The apparent disagreement between the dramatic devi- 
ation of the ratio from unity seen in the EMC data and 
the small nuclear effects predicted by theoretical calcula- 
tions has triggered a series of further measurements and 
theoretical investigations, for reviews, see 041] • 

The emerging picture of nuclear modifications of the 
nucleus to deuteron cross section ratio, a A /a D , has the 
pattern presented in Fig. [T] For small values of Bjorken 
x, x < 0.05 — 0.1, the ratio is noticeably suppressed — 
the suppression increases with an increase of the atomic 
number A and a decrease of x — which is called nuclear 
shadowing. For 0.1 < x < 0.3, the ratio is enhanced; 
the effect is small (of the order of a few percent) and 
does not reveal an obvious A dependence. In the interval 
0.3 < x < 0.8, the ratio is suppressed and this suppres- 
sion is called the EMC effect. Finally, for x > 0.8 the ra- 
tio dramatically grows above unity which is explained by 
the effect of the nucleon motion inside nuclei (Fermi mo- 
tion). Various models describe the experimental a A /<j D 
cross section ratios for certain ranges of Bjorken x, but 
there is no comprehensive understanding of the entire 
pattern of the nuclear modifications described above. In 
particular, there is no unique and generally accepted the- 
ory to explain the nature of the antishadowing and EMC 
effects. 

In this paper we focus on the enhancement (antishad- 
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FIG. 1: The pattern of nuclear modifications of the a A /a D 
cross section ratio as a function of Bjorken x for 56 Fe and 
64 Cu. The data are from BCDMS 01 (open squares), SLAC 
E139 0| (filled circles) and EMC [1] (stars). For all data 
sets statistical and systematic errors have been combined in 
quadrature. 



owing) of the a ja D cross section ratios in the 0.1 < x < 
0.3 region. The deviation of a A /a D from unity in the an- 
tishadowing region is of the order of a few percent 
(see Fig. [J) . Given that most measurements quote nor- 
malization uncertainties on the order of 1-2% (usually 
due to target thickness or luminosity), it is difficult to 
quantify the absolute size of the antishadowing effect pre- 
cisely, and comparisons between experiments are some- 
what complicated. In addition, systematic uncertainties 
due to radiative corrections are highly non-trivial in this 
region of x, and are sometimes hard to determine accu- 
rately. An example of the difficulty involved in achieving 
very precise measurements in the antishadowing region 
can be found in the SLAC E139 results. The preliminary 
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results for the Fe/D ratio were essentially consistent with 
unity for the region 0.1 < x < 0.3 @- However, the final 
E139 analysis yielded results in the antishadowing region 
more consistent with, e.g., the EMC and BCDMS exper- 
iments, showing a small enhancement of « 3% on aver- 
age • Despite the difficulties inherent in antishadowing 
measurements, the results from various experiments are 
remarkably consistent within their experimental uncer- 
tainties. In addition, the small enhancement seen by the 
EMC, BCDMS, and SLAC E139 experiments for cop- 
per and iron targets has also been seen in lighter targets 
(Ca/DjN/D, C/D, He/D) by the NMC 10} and HER- 
MES [lj experiments. 

The antishadowing effect has rather intriguing fea- 
tures. Unlike the shadowing effect, antishadowing 
showed little or no sensitivity to the mass number A 
within experimental uncertainties, for example, in the 
SLAC E139 and NMC data 0- While antishadow- 
ing is observed in nuclear DIS, the cross section enhance- 
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ment is not seen in nuclear Drell- Yan rates 12 1 and total 
neutrino-nucleus cross sections for x > 0.1 13 1. 

In the leading twist formalism, the small enhancement 
of a A J a D in the antishadowing region translates into an 
enhancement of the valence quark and possibly gluon dis- 
tributions in nuclei in this region |13T[l7| . However, the 
pattern and especially the magnitude of nuclear modifi- 
cations of the gluon distribution in nuclei are very poorly 
constrained by present data. 

The aim of this paper is to examine the impact of the 
nuclear dependence of R = (tl/&t, i.e., the ratio of the 
longitudinal to transverse photoabsorption cross sections, 
on the extraction of the nucleus to deuteron structure 
function ratios, F A /F® and Ff/Ff, from a A /a D data. 
In particular, we demonstrate that in the presence of a 
small but non-zero difference between R for nuclei and 
the nucleon, the nuclear enhancement in the ratio of the 
transverse structure functions F A /F® becomes signifi- 
cantly reduced (or even disappears in some cases), indi- 
cating that antishadowing is dominated by the longitudi- 
nal contribution. In addition, we analyze how the nuclear 
dependence of R affects the nuclear gluon distribtion and 
emphasize the importance of measurements of R in the 
DIS kinematics as a direct probe of the gluon distribution 
in nuclei. 



2. NUCLEAR DEPENDENCE OF R AND THE 
RATIO OF NUCLEUS AND DEUTERON 
STRUCTURE FUNCTIONS 

2.1. Longitudinal contribution to the inclusive 
cross section 

In the one-photon exchange approximation, the spin- 
independent cross section for inclusive electron scattering 
can be expressed as 
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where ot (<xl) is the cross section for photoabsorption 
of purely transversely (longitudinally) polarized photons; 
R = <7l/&t] r is the transverse virtual photon flux; 
e is the virtual photon polarization parameter. In the 
laboratory frame, the negative four-momentum squared 
(virtuality) of the exchanged photon is — q 2 = Q 2 = 
4EE'sm 2 (8/2) and the Bjorken x is x = Q 2 /[2M(E - 
E')], where E (£") is the energy of the incident (scat- 
tered) electron, 9 is the scattering angle, and M is the 
nucleon mass. The flux of transverse virtual photons can 
be expressed as T = aE'(W 2 ~ M 2 )/{4tt 2 Q 2 ME(1 - e)], 
where a is the fine structure constant and W is the invari- 
ant energy of the virtual photon-proton system. Finally, 
the virtual photon polarization parameter is: 
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where v = E — E'; y = vjE. Note that in the second line 
of Eq. ([2]) we expressed e in a Lorentz invariant form. 

In terms of the structure functions Fi(x,Q 2 ) and 
F2{x, Q 2 ) in the DIS region, the double differential cross 
section can be written as 
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A comparison of Eqs. ([T]) and §3§ shows that Fi(x, Q 2 ) 
is purely transverse, while 

AM 2 x 2 

F L {x,Q 2 ) = (l + —^)F 2 (x,Q 2 )-2xF 1 (x,Q 2 ) (4) 

is purely longitudinal. Note that F 2 (x, Q 2 ) is a mixture 
of both the longitudinal and transverse contributions. 
Thus, the ratio R is 



R 



ol = F L (x,Q 2 ) 
a T 2xF 1 (x,Q 2 ) 



(5) 



The nucleon structure function F 2 (x,Q 2 ) is propor- 
tional to the d 2 a/(dildE') differential cross section in 
the e — > 1 limit; it has been measured with high preci- 
sion in various x and Q 2 bins. The longitudinal structure 
function Fl(x, Q 2 ), in contrast, is not measured as well as 
F 2 (x, Q 2 )'- the data is sparse and imprecise for the proton 
and is even more limited for nuclei. It is an experimental 
challenge to separate F 2 (x,Q 2 ) and Fl{x,Q 2 ) which is 
usually done using the method of Rosenbluth separation, 
i.e., by measuring the cross section at different energies 
(at fixed x and Q 2 ) to allow for a variation of e. 
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In this paper we shall use the parameterization of R 
for the nucleon, R N , that was obtained as the result of 
the global analysis of the SLAC hydrogen and deuterium 
data . The same analysis also showed that R D = R H 
to high accuracy, where R D (R H ) refers to deuterium 
(hydrogen) . An example of the values of R N in the kine- 
matics used in this paper is presented in the middle panel 
of Fig. [21 Note that the more recent analysis of the SLAC 
E143 collaboration J19] reported results for R N consistent 
with those of Ref. [3- 



2.2. Hints of nontrivial nuclear dependence of R 

Experimentally measured cross section ratios contain 
both transverse and longitudinal contributions of the 
structure functions. In terms of the structure function 
F2(x,Q 2 ), one can write the ratio of the nucleus to 
deuteron photoabsorption cross sections as 
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where the superscript A refers to the nucleus and the 
superscript D refers to the deuteron; AR = R A — R D . In 
the second line of Eq. ^ we used the expansion in terms 
of the small parameter AR and kept first two terms of 
the expansion. 

Alternatively one can express the cross sections a A and 
a D in terms of the structure function F\(x, Q 2 ) and ob- 
tain: 
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The cross section ratio a A ja D can be identified with the 
structure function ratio F A / F® or F A /FF only with 
the assumption of the trivial nuclear dependence of R = 
cl/ct, he., R A — R D , or in certain kinematic limits. In 
particular, a A /<r D = F A / F® at e = 1 and a A ja D = 
F A /Ff at e = 0. 

Figure [2] presents the kinematic coverage in Q 2 and x 
and the corresponding values of e and R N of the data 
points considered in this paper. On the one hand, the 
BCDMS [|, EMC H and NMC [H| data are mostly 
taken with e close to unity (see the upper panel of Fig. [5]), 
which implies that the cross section ratios are close to the 
F2 structure function ratios, even if AR = R A — R D ^ 0. 
On the other hand, the SLAC data 0, Hq| (see the up- 
per panel of Fig. [2]) corresponds to the kinematics where 
e « 0.5 and, hence, F A jF® will deviate from a A /a D if 
AR ^ 0. For all these experiments t / and, hence, 
the extraction of the transverse structure function ratios 
F A I F-f depends explicitly on the assumption adopted 
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FIG. 2: The kinematic coverage in Q and x and the corre- 
sponding values of e and R of the data points considered in 
this paper. 



for AR. Below we summarize what is known about it 
from the world data. 

At small Q 2 , R might be different between deuterium 
and hydrogen [2ll|. though it seems to be identical at 
large Q 2 [22|, |23j. In particular, there are some hints in 
both JLab E99-118 and SLAC data [H that R D 
is smaller than R H for Q 2 < 1.5 GcV 2 , with a global 
average of R D - R H = -0.054 ± 0.029. 

Turning to heavier nuclei, the SLAC E140 data [Io| 
suggest some nuclear dependence of R at x — 0.2, which 
seems to have a nontrivial Q 2 dependence: R Fe — R D can 
be positive at Q 2 = 2.5 GeV 2 and negative at Q 2 = 1.5 
and 1 GeV 2 : 



R Fe 
R Fe 
R Fe 



P=2.5 



R D \ 
-R D |q2 = 1. 5 

R 



0.144 ± 0.079(stat.) ± 0.027(syst.); 
-0.124 ± 0.051(stat.) ± 0.023(syst.): 
-0.086 ± 0.057(stat.) ± 0.022(syst.). 
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A word of caution is in order here. Coulomb corrections 
may be non-negligible in DIS at SLAC and JLab kine- 
matics, especially at large x. These corrections tend to 
reduce R for nuclear targets [24} ■ 

The nuclear dependence of R at Q 2 of the order of 
a few GeV 2 and lower was also measured by the HER- 
MES collaboration by fitting the a A /a D cross sec- 
tion as a function of the virtual photon polarization e 
over a typical range of 0.4 < e < 0.7. Overall no signifi- 
cant nuclear dependence of R for 14 N and 3 He targets for 
x > 0.06 has been observed (the data in the x < 0.06 re- 
gion is affected by the correlated background and should 
be neglected). However, since this was a single-energy 
measurement with correlated values of e and Q 2 , the ex- 
traction of R A /R D was done in a mo del- dependent way. 

At larger values of Q 2 , the NMC experiment [25| ob- 
tained R Ca - R c = 0.027 ± 0.026(stat.) ± 0.020(syst.) 
at (Q 2 ) — 4 GeV 2 and concluded that AR is compati- 
ble with zero. However, a hint of the nontrivial nuclear 
dependence of R can be still seen in data. The preci- 
sion Sn/C data from NMC dl show that R Sn - R c = 
0.040 ± 0.021(stat.) ± 0.026(syst.) at a mean Q 2 of 10 
GeV 2 . This value of AR = R A - R D corresponds to 
AR/R N = 0.22 + 1.20, i.e., 22 - 120 % relative devia- 
tion for different values of x in the considered kinematics, 
where R N is given by the parameterization of Ref. [l8[ 
presented in the middle panel of Fig. [5J Note that the ex- 
traction of AR in this experiment was based on a method 
closely related to Rosenbluth separation taking advan- 
tage of three different incident muon energies (120, 200 
and 280 GeV). 

For convenience, we present in Table U a brief overview 
(covered kinematics in Bjorken x and Q 2 and energy set- 
tings) of the discussed measurements of the nuclear de- 
pendence of R (involving nuclei heavier than deuterium) . 
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FIG. 3: R Sn -R c as a function of x. Full squares are results of 
the NMC measurement with the statistical and systematic er- 
rors added in quadrature [26j] ; the long-dash and dotted curves 
correspond to R Sn - R c = 0.04 and R Sn - R c = 0.3R N , re- 
spectively; the curves labeled "EPS09" , "HNK07" and "nDS" 
correspond to predictions using different nuclear parton dis- 



tributions. Also shown are the NMC result for R' 



Ca 



R L 



(triangle), the SLAC result for R Au - R' 



gleh and SLAC E140 results for 
x [2(j] (open circles). 



R" 



e M (inverse trian- 
R D as a function of 



Experiment and 
observables 


Kinematics 


Beam energy 


SLAC E140 [20] 

j^Fe,Au j^D 


0.2 < x < 0.5 
1 < Q 2 < 10 GeV 2 


3.7 < E < 15 GeV, 
up to five energies 


NMC (1992) [25] 

R Ca _ R C 


0.0085 < x < 0.15 
1 < Q 2 < 15 GeV 2 


E = 90 and 200 GeV, 
two energies 


HERMES [11] 

R *He,Ne /R D 


0.01 < x < 0.8 
0.2 < Q 2 < 3 GeV 2 


E = 27.5 GeV 
single energy 


NMC (1996) [26J 
R Sn - R c 


0.0125 < x < 0.45 
3.3 < Q 2 < 35 GeV 2 


E = 120, 200, 280 
GeV, three energies 



TABLE I: An overview of the measurements of the nuclear 
dependence of R discussed in this paper. 

The results of the NMC measurement of R s " —R c as a 
function of Bjorken x (26| are presented as full squares in 
Fig. [3J For completeness, we also show the NMC result 
for the average R Ca 



R 



c 



25[ as a trian gle, the SLAC 
E140 result for the average R Au - R Fe [27] as an in- 
verse triangle, and the SLAC E140 results for R Fe - R D 
as a function of x : 20] as open circles. (For the lat- 
ter, we showed only the data points for the 6% radiation 



length iron target and shifted them in x for better visi- 
bility.) The long-dash and dotted curves correspond to 
R Sn _ R c = o 04 and R sn _ R c = q^R n , respectively. 

As one can see from the figure, both curves provide a 
good description of the measured values of R Sn — R c . Fi- 
nally, the curves labeled "EPS09" , "HNK07" and "nDS" 
correspond to the direct calculation of R Sn — R using 
different parameterizations of leading twist nuclear par- 
ton distributions (PDFs), see the discussion in Sect. H] 
Note that we have singled out the NMC Sn/C data [26[ 
because the extraction of R A — R D was done using a 
method closely related to the Rosenbluth separation and 
because the covered kinematics (the values of x, Q 2 and 
e) broadly overlaps with that of the BCDMS, EMC, and 
NMC data on <r A ja D that we analyze in this paper. 

In summary, as a global average, while R seems to 
show little nuclear dependence within relatively large ex- 
perimental uncertainties, there exist hints of nontrivial 
nuclear dependence of R. In particular, AR = R A — R D 
may be statistically different from zero in some kinemat- 
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Impact of nuclear dependence of R on nucleus 
to deuteron structure function ratios 



2.3. 



As we explained in Sect. 12.21 if there is a nontriv- 
ial nuclear dependence of R, the a A /a D cross section 
ratio is not equal to the F A / Ff or F A / F® structure 
ratios. In particular, a positive R A — R D will lead to 
Ff/Ff < a A /a D < F A /F®. Since the nuclear depen- 
dence of R has not as yet been systematically measured, 
we shall test two assumptions for AR that are motivated 
purely by the NMC Sn/C data [26[, which has kinematic 
coverage similar to that of the BCDMS, EMC and NMC 
measurements. In our analysis below we assume that: 

1) (Absolute) AR = R A - R D = 0.04. This is based 
on the NMC measurement of R Sn — R c at an average 
(Q 2 ) = 10 GeV 2 . 

2) (Relative) (R A -R D )/R N = 30%, which is possible in 
view of the fact that the NMC Sn/C data allows for the 
22 - 120% relative deviation of AR/R N . 

Note that we effectively assumed that R A — R D « 
R Sn — R c which corresponds to the lower limit for AR. 

The impact of our assumptions for AR on selected nu- 
clear DIS data is presented in Figs. [4] and [5j we truncated 
the used data sets by neglecting low x and high x data 
points and focusing on the antishadowing region. 

The BCDMS Fe/D [6], EMC Cu/D [8] and NMC 
Ca/D [Io( data presented in Fig. [4] correspond to e close 
to unity. Therefore, regardless of the assumption for AR, 
one expects that F A jF® w <t a /<t d with a very good ac- 
curacy. On the other hand, F A / 'Ff is clearly smaller 
than g a ja D . Thus, the few percent enhancement of 
a A /a D in the antishadowing region may be reduced or 
removed altogether for the ratio of the transverse struc- 
ture functions F A /Ff. 

For the SLAC E139 Q and E140 [H Fe/D data pre- 
sented in Fig. [5j the values of Q 2 are rather small (see 
the lower panel in Fig. [2]) and our assumptions for the 
nuclear dependence of R motivated by the NMC Sn/C 
measurement at higher Q 2 require a significant extrapo- 
lation in Q 2 . However, for the lack of better input, in our 
analysis of the SLAC data we adopt the same assump- 
tions for AR as those used above. Since the values of e for 
these two data sets are not close to unity (see the upper 
panel in Fig. [5]) , AR > leads to noticeable differences 
between the ratio of the structure functions and the ratio 
of the cross sections according to the trend described by 
Eqs. © and ©: F A /Ff < a A /a D < F A /F®. 

In summary, the assumed nontrivial nuclear depen- 
dence of R leads to a decrease or to a complete disap- 
pearance (in some case) of enhancement of the F A /Ff 
structure function ratio in the 0.1 < x < 0.3 region. If 
confirmed by future experiments, this observation would 
indicate that the effect of antishadowing in a A is pre- 
dominantly due to the contribution of the longitudinal 
structure function F A , instead of F A as implicitly as- 
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FIG. 4: The impact of the nontrivial nuclear dependence of 
R on the structure function ratios around the antishadow- 
ing region for BCDMS Fe/D @|, EMC Cu/D @ and NMC 
Ca/D [lCj data. The values of e are close to unity. 



sumed in most phenomenological analyses. 



EXPERIMENTAL LIMITS ON 



DETERMINING 7?' 



R L 



Thus far we have examined the impact of a nuclear de- 
pendence of R on the extraction of the nuclear dependent 
structure function ratios F A /Ff and F A /F® from cross 
section ratios. The logical question then becomes: 'What 
is the limit on the experimental precision for R A — 
In this section we shall explore this question within the 
context of the precision likely to be avialable for dedi- 
cated L/T separation measurements over the next decade 
or two. For guidance we shall refer to the highest pre- 
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FIG. 5: The impact of the nontrivial nuclear dependence of 
R on the structure function ratios around the antishadowing 
region for SLAC E139 and E140 HJ Fe/D data. 




FIG. 6: Simulated <r A /a D as a function of e (data points 
with error bars) and a linear fit (solid lines), see the text for 
details. 



cision experiments performed at SLAC @, Hoj, HH and 
Jefferson Lab [U [28[. These experiments have shown 
that reducing the a a/vd cross section ratio uncertainties, 
point-to-point in e, below 1% is experimentally challeng- 
ing, yet obtainable. For instance, the point-to-point un- 
certainties from Jefferson Lab experiment E94110 [28| on 
cryogenic hydrogen have been estimated at about 1.5%, 
which was found to be consistent with the width of the 
distribution of residuals determined from the linear fits. 

To measure cross sections at a range of e values for 
fixed x and Q, both the SLAC and JLab inclusive L/T 
separation experiments utilized a range of beam energies 
in conjunction with well studied spectrometer systems, 
which were able to be rotated to different angles and ad- 
justed to accept varying ranges of momenta. Some of 
the largest contributions to the estimated systematic un- 
certainties stem from either time dependent systematics, 
such as current calibrations or detector efficiency varia- 
tions, or from the uncertainties in the kinematics at each 
beam energy and spectrometer setting. However, these 
systematics largely cancel in the cross section ratios, in 
which the electron yields on each target are taken at the 
same kinematic settings and close in time. 

If, for example, a 3% anti-shadowing effect in F A 
were entirely due to a longitudinal enhancement, with 
F A /F® = 1, then this would be reflected in a 3% slope 
in the cross section ratio versus e' = e/(l + eR D ), cor- 
responding to AR — R A — R D w 0.03. For the current 
study we assume the following: 



• The total systematic point-to-point uncertainty (in 
e) on the measured a A ja D ratios is 0.5%. 

• There is no e dependent systematic uncertainty. 

• Six cross section ratio measurements at equally 
spaced e' = e/(l + eR D ) values in the range (0.05, 
0.95), corresponding to six unique beam energies. 

Under the assumptions above, the cross section ratios 
were selected at each e' by random sampling from a Gaus- 
sian distribution assuming a 3% slope on a A /a D versus 
e' and a Gaussian width of 0.5%. Six sample L/T sepa- 
rations generated by this procedure are shown in Fig. [5] 
After performing a linear fit, the uncertainty on the mea- 
sured slope was found to be 0.67%, corresponding to a 1-a 
(3-cr) uncertainty on R A - R D of less than 0.007 (0.021). 
For the case considered of 0.5% ratio uncertainties, one 
could determine at 1-a whether a 3% antishadowing ef- 
fect is due mainly to F A to « 20%. 

We note that this uncertainty on the extracted R A — 
R D scales with the uncertainties on the cross section ra- 
tios such that a further reduction in the latter to 0.25% 
would reduce the uncertainty on AR by half. However, 
we have thus far ignored any possible e dependent sys- 
tematic uncertainties, such as those possibly arising from 
Coulomb and radiative corrections. For this reason, this 
is likely an optimistic scenario. 
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4. NUCLEAR DEPENDENCE OF R AND ITS 
ROLE IN ANTISHADOWING OF THE GLUON 
DISTRIBUTION IN NUCLEI 



We demonstrated in Sect. 12.31 that the assumption 
of the nontrivial nuclear dependence of R, i.e., R A — 
R D > 0, whose magnitude and sign are motivated 
by the NMC Sn/C data [Hj], leads to a difference be- 
tween the cross section and structure function ratios: 
Ff/Ff < <r A /a D < F A /Ff. Moreover, the reduction 
of the F A I F® ratio is quite sizable: the enhancement in 
the 0.1 < x < 0.3 region visible in the cross section ra- 
tios is significantly decreased (or even disappears) for the 
F A /F® ratios, which indicates that antishadowing pre- 
dominantly resides in the longitudinal structure function 
F A . This conclusion is rather general; in particular, it 
does not rely on the twist expansion and the underlying 
partonic structure. 

In the framework of the leading twist formalism, global 
QCD fits to the available data [l3|-[lj| show that the 
small enhancement of a A ja D in the antishadowing re- 
gion translates into an enhancement of the valence quark 
and possibly gluon distributions in nuclei compared to 
those in the free proton. One should emphasize that 
these analyses assumed no nuclear dependence of i?, i.e., 
R = 0. The pattern and magnitude of nuclear modifi- 
cations of the nuclear gluon distribution gA (x) is known 
with very large uncertainty because gA{%) is mostly de- 
termined indirectly from scaling violations of the nuclear 
structure function F A measured in a limited kinematics. 
This is illustrated in Fig. [7J where we present the ratio 
of leading order gluon distributions in 40 Ca to that in 
the free proton, gA (x) / g^ (x) , as a function of x at fixed 
Q 2 = 3 GeV 2 . In the figure, the solid curve is the result 
of the EPS09 fit [lj]; the dotted curve is the result of the 
HKN07 fit 1151 : the dot-dashed curve is the nDS parame- 
terizationjlj] , whose results are quantitatively similar to 
those of [l^For the EPS09 and HKN07 fits, we showed 
only the central values; the theoretical uncertainty on 
these predictions is quite large essentially in the entire 
range of x. 

As one can see from Fig. [JJ different groups predicts 
wildly different gA{x)/giq{x) (with large uncertainties). 
Since the amounts of nuclear shadowing and antishad- 
owing are correlated through the momentum sum rule, 
large antishadowing corresponds to significant shadow- 
ing in the EPS09 fit jiff]; very small antishadowing corre- 
sponds to negli gibly small shadowing in the nDS fit [Hj]); 
the HKN07 fit [15[ suggests yet another scenario where 
large gluon antishadowing is concentrated at large x. 

Given the present uncertainty in gA (x), it is impor- 
tant and instructive to confront the NMC measurement 
of R Sn — R c [26| with direct calculations of this quantity 
in the framework of leading twist nuclear parton distri- 
butions. This is presented in Fig. [3] where the curves 
labeled "EPS09", "HKN07" and "nDS" correspond to 
the direct calculation of R — R in the kinematics of 
the NMC measurement [26] using the respective leading 



x 



< 

CD 




FIG. 7: gA(x)/gN(x) for 40 Ca as a function of x at fixed 
Q 2 = 3 GeV 2 as obtained from global QCD fits. The solid 
(dotted, dot-dashed) curve is the result of the EPS09 QiJ 
(HKN07 fit nDS 0) fit. 



order parton distributions in nuclei. One can readily see 
from the figure that while for small x, x < 0.05, the 
leading twist description is consistent with our assump- 
tions for Ai? 7^ and the NMC data, in the antishad- 
owing region 0.1 < x < 0.3 and also for larger x the 
leading twist approach predicts a negligibly small Ai? 
in contrast with our assumptions and only marginally 
agrees with the data due to the large experimental un- 
certainty. Note that the leading twist calculations pre- 
sented in Fig. [3] have rather small theoretical uncertain- 
ties stemming mostly from the uncertainty in the gluon 
distributions. 

There are several reasons for the negligibly small value 
of R Sn - R c for x > 0.05 at the NMC energies predicted 
in the leading twist framework. First and most impor- 
tantly, the assumed shapes of the parameterizations of 
quark and gluon distribution in nuclei [l4T - fl7| are such 
that nuclear PDFs and the ratio R show only a weak 
nuclear dependence around x — 0.1 (see Figs. [7J and [5]). 
For instance, while the EPS09 analysis [l6[ used the data 
on Q 2 dependence of F^ n (x,Q 2 )/Ff(x,Q 2 ) [f|, it did 
not include the R — R data in the fit. Hence, re- 
sulting nuclear PDFs were not constrained to reproduce 
the experimental values of R Sn — R c which, as a re- 
sult, leads to R Sn - R c w for x > 0.1. Second, while 
R Sn /R N and R c / R N separately reveal quite sizable de- 
viations from unity (compare to R Ca /R N presented in 
the upper panel of Fig. [5]), nuclear effects mostly cancel 
in the R Sn — R c difference. In general, while it is natural 
to expect Ai? ^ because the pattern of nuclear modifi- 
cations of quark and gluon distributions is different, with 
the currently assumed shapes of nuclear parton distribu- 
tions it is not easy to generate sizable Ai? for x > 0.1 
and large Q 2 because i? itself is very small there. Third, 
in the NMC kinematics the values of Bjorken x > 0.1 
correspond to Q 2 > 10 GeV 2 . At such large values of 
Q 2 , nuclear modifications of parton distributions gradu- 
ally become less pronounced. Note also that it is unlikely 
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that higher twist (twist-four) effects can generate sizable 
AR because it would require unrealistically large higher 
twist effects 0. 

While the available data on the nuclear dependence of 
R is not able to constrain the nuclear gluon distribution 
in the 0.1 < x < 0.3 region, a better chance of measuring 
gluon antishadowing would be offered by measurements 
of R with nuclear targets and the deuteron (proton) and 
at not too high Q 2 . Note this is essentially equivalent to 
measuring the longitudinal structure functions Fl(x, Q 2 ) 
for nuclei and the deuteron (proton). Such measurements 
can be carried out at Jefferson Lab at 12 GeV at low-to- 
intermediate Q 2 [iH and at a future Electron-Ion Collider 
(EIC) at intermediate-to-high Q 2 (HSH. In the latter 
case, the measurement of F A (x, Q 2 ) and the longitudinal 
nuclear structure function F A (x,Q 2 ) (taking advantage 
of variable energies) with the subsequent extraction of 
qa (x) in a wide kinematic range is already an important 
part of the planned physics program. 

An example of expected nuclear effects is presented in 
Fig.[8]which shows predictions for the ratio of the nuclear 
to nucleon ratios R A /R N (upper panel) and longitudinal 
structure functions F A / F^ (lower panel) as a function 
of x at Q 2 = 3 GeV 2 for 40 Ca. Different curves cor- 
respond to different parameterizations of nuclear PDFs 
(see Fig. [7]). A comparison of Figs. [7J and [5] shows that 
different assumptions about the shape of the gluon (and 
quark) distributions in nuclei lead to different shapes of 
R A /R N and F A /F^ . To point out just one feature, an 
observation of sizable R A /R N > 1 (enhanced F A /F^ 
compared to Ff/Ff) and F A /F^ > 1 in the antishad- 
owing region 0.1 < x < 0.3 would unambiguously signal 
the presence of a significant antishadowing for the gluon 
distribution in nuclei. (The gluon distribution enters the 
longitudinal structure function Fl(x, Q 2 ) at the same or- 
der as the quark distributions; at the same time, the 
gluon distribution enters the transverse structure func- 
tion F\(x, Q 2 ) with the weight (coefficient function) that 
is smaller than that for ^(cc, Q 2 ) [HJ .) The converse is 
also true: an absence of nuclear enhancement of R A /R N 
and F A /Fj^ in the interval 0.1 < x < 0.3 would trans- 
late into the absence of antishadowing for gluons in this 
region. 



5. CONCLUSIONS 

In this paper we studied the influence of the nontriv- 
ial nuclear dependence of R = <jl/^t on the extrac- 
tion of the F A jF^ and F A / F® structure function ra- 
tios from the data on the a A /a D cross section ratios. 
Guided by indications of the nuclear dependence of R 
from the world data and, in particular, by the NMC 
measurement that showed that R Sn — R c = 0.040 ± 
0.021(stat.) ± 0.026(syst.) at (Q 2 ) = 10 GeV 2 [Hj], we 
tested two assumptions for AR = R A - R D : AR = 0.04 
and AR/R N — 0.3, where R N corresponds to the free 
proton [18j. With these assumptions, we examined se- 
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FIG. 8: R A /R N (upper panel) and F A /F^ (lower panel) 
qa{x) I 'gN(x) for 40 Ca as functions of x at Q 2 = 3 GeV 2 . The 
solid (dotted, dot-dashed) curve is the result of the EPS09 [l6l ] 
(HKN07 fit 111, nDS [Hf) fit. 



lected sets of EMC, BCDMS, NMC and SLAC data on 
a A /a D and extracted the F^/F? and Ff/Ff ratios. 
We find that for the EMC, BCDMS and NMC data, 
F^/Ff « a A /a D , while Ff/Ff < a A /a D . For the 
SLAC data, we found that F A /Ff < a A ja D < F A /F°. 
In particular, we observed that the nuclear enhance- 
ment (antishadowing) in the interval 0.1 < x < 0.3 be- 
comes significantly reduced (or even disappears in some 
cases) for the ratio of the transverse structure functions 
F A I F® . The latter observation indicates that antishad- 
owing may in fact be dominated by the longitudinal con- 
tribution rather than by the transverse one (i.e., anti- 
shadowing is dominated by gluons rather than by quarks) 
as implicitly assumed by current phenomenological anal- 
yses and global nuclear parton distribution fits. 

We also examined experimental limits on determining 
R A — R D from measurements of the e' — (-/{! + eR D ) de- 
pendence of a A I a D . Making a plausible assumption that 
a A /a D has a 3% slope in e' and can be measured with a 
0.5% uncertainty over a broad range of e', we found that 
AR can be extracted with 0.67% uncertanity. Therefore, 
one could determine whether a 3% antishadowing effect 
is mainly due to F A to approximately 20% accuracy. 

In the leading twist framework, the magnitude of nu- 
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clear enhancement of R A and the longitudinal struc- 
ture function F^(x,Q 2 ) (these quantities directly probe 
the nuclear gluon distribution gA{x)) is directly corre- 
lated with the size and shape of antishadowing for gA(x). 
While at the moment gA (x) is rather poorly constrained 
by QCD fits to available data, a dedicated high-precision 
measurement of the nuclear dependence of R (the longi- 
tudinal nuclear structure function Ff{x,Q 2 )) at Jeffer- 
son Lab and an EIC has the potential to unambiguously 
constrain gA (x) in the antishadowing region and beyond 
(An EIC will also be able to constrain gA (x) deep in 
the shadowing region of small x.) Through the parton 
momentum sum rule, this knowledge will have a deep 
impact on <?a(x) in the entire range of x. In particular, 
it should dramatically help to constrain gA{x) in the nu- 



clear shadowing region, < x < 0.05, where gA{x) 

plays an essential role in phenomenology of high-energy 
hard processes with nuclei, for a review, see [33j ]. 
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